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Abstract. We investigate the translocation of a stiff polymer consisting of M 
monomers through a nanopore in a membrane, in the presence of binding particles 
(chaperones) that bind onto the polymer, and partially prevent backsliding of the 
polymer through the pore. The process is characterized by the rates: k for the polymer 
to make a diffusive jump through the pore, q for unbinding of a chaperone, and the rate 
qn for binding (with a binding strength k) ; except for the case of no binding k — Q the 
presence of the chaperones give rise to an effective force that drives the translocation 
process. In more detail, we develop a dynamical description of the process in terms 
of a (2+1) variable master equation for the probability of having m monomers on the 
target side of the membrane with n bound chaperones at time t. Emphasis is put on 
the calculation of the mean first passage time j as a function of total chain length 
M . The transfer coefficients in the master equation are determined through detailed 
balance, and depend on the relative chaperone size A, binding strength k, as well as 
the two rate constants k and q. The ratio 7 = q/k between the two rates determines, 
together with n and A, three limiting cases, for which analytic results are derived: (i) 
For the case of slow binding (7K 0), the motion is purely diffusive, and 7 ~ for 
large M; (ii) for fast binding (7K 00) but slow unbinding (7-^0), the motion is, for 
small chaperones A = 1, ratchet-like, and j ~ Af; (iii) for the case of fast binding and 
unbinding dynamics (7 ^ 00 and 7K 00), we perform the adiabatic elimination of 
the fast variable n, and find that for a very long polymer j ~ M, but with a smaller 
prefactor than for ratchet-like dynamics. We solve the general case numerically as a 
function of the dimcnsionless parameters A, k and 7, and compare to the three limiting 
cases. 
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1. Introduction 

The passage of a biopolymer (DNA, RNA, proteins) tlirougli a nanopore embedded in 
a membrane is one of the most crucial processes in biology, examples including: the 
translocation of proteins through the endoplasmatic reticulum, the passage of RNA 
through the nucleus pore membrane, the viral injection of DNA into a host, as well as 
DNA plasmid transport from cell to cell through cell walls pj. Moreover, biotechnology 
applications of translocation processes such as rapid DNA sequencing [21IS1, secondary 
structure determination of RNA f^, as well as analyte detection and nanosensing [30111 
have been discussed. Goals such as drug delivery are among the ultimate research 
directions connected to translocation. 

Experimentally, biopolymer translocation can be probed on the single molecule 
level. In such experiments one commonly employs an a-hemolysin protein pore 
introduced into a planar lipid bilayer. But, of late, nanopores can also be created in a 
well-defined manner by etching techniques in silicon oxide |Hj. In vitro, translocation 
experiments involve either a constant driving voltage maintained during the entire 
measurement [HlEl; or a high applied voltage for threading the head of the biopolymer 
through the pore, followed by an off- voltage and a low probe voltage [Ej. Once a 
biopolymer enters the pore, it blocks the electrical current of ions through the pore, as 
monitored by patch-clamp techniques [Oj. From such measurements one may directly 
extract, for instance, the mean translocation time j and its dependence on applied 
voltage and on other physical parameters of the system. Note that cross-membrane 
potentials driving translocation also exist in nature ^J. 

Another in vivo mechanism is the driving by binding proteins, so-called chaperones, 
that appears to be particularly common for protein translocation [TT | IT ^ IT! ^ ITi t ITT f IH?] . 
but is also implicated in the DNA transport through membranes [T^ I2Uj. The driving of 
the translocation process by binding proteins constitutes a fine example of how directed 
motion can emerge from the coupling of two random processes. Schematically, this is 
illustrated in figure E binding proteins on the trans side (B) of the membrane bind 
nonspecifically to the already translocated portion of the biopolymer. One expects 
three limiting dynamical regimes: (i) for slow binding (in reference referred to as 
the diffusive regime), the relative diffusive motion of the polymer through the pore is 
so fast that the chaperones do not have time to bind; (ii) for fast binding but slow 
unbinding (the irreversible binding regime), the motion is ratchet-like: The chaperones 
bind (infinitely) fast and while being attached to the polymer, they prevent backsliding 
through the pore, and therefore (partially) rectify the diffusive process, (iii) for the 
case of fast binding and unbinding dynamics (the reversible binding regime), one may 
assume that the (fast) chaperones exert an effective force originating from the chaperone 
chemical potential difference between the trans and cis sides |22l 1231 123 12^1 12E] • Given 
the rather short biopolymer segments used in typical experimental setups as well 
as short translocating biopolymers in vivo jl], in combination with the rather large 
binding interface between chaperones and driven biopolymer, non-negligible finite size 




Figure 1. Schematic of the translocation process: Binding proteins (chaperones) , of 
volume concentration cq, attempt to bind to the translocating chain. Once bound, 
chaperones (partially) rectify the passage of the biopolymer through the nanopore. 
The strength of the binding between the polymer and the chaperones is characterized 
by the dimensionless parameter k, = cqK^'^, where K""^ is the chaperone binding 
constant, see text. When bound onto the polymer a single chaperone occupies a 
length Act, where a is the linear monomer size (size of a base, base-pair, or aminoacid). 
The chaperones span several monomeric unit lengths, typically A ^ 10 — 20. In the 
absence of chaperones the polymer diffuses through the pore with a diffusive rate k 
(the bare diffusion constant for the chain is thus Do = c^fc). The (un)binding process 
is characterized by an unbinding rate q for a single chaperone, leading to a competition 
of two random processes depending on the ratio ^ — q/k. 



effects occur, that require a modeling beyond the effective chemical potential. Such 
finite size effects were studied in detail in reference [21], finding distinct oscillations in 
the effective driving force. Although the different dynamical regimes above have been 
investigated in detail, so far no model has been able to account for all regimes and 
the switching in between them, as tuned by changing the physical parameters, that is 
possible in in vitro experiments. A first step towards such a unification was taken in 
reference jTTj where Brownian dynamics simulation were performed, and compared to 
the results obtained from a (2-|-l)-variable Fokker-Planck equation. 

We here extend the work in reference ^1] and propose a (2-|-l)-variable master 
equation formalism for describing the full coupled dynamics of chaperone (un) binding 
and polymer diffusion during the translocation. In contrast to the continuum 
description in reference our discrete approach can explicitly account for the ratchet 
mechanism, and include chaperones which are larger than the size of a monomer. Our 
approach resembles recent models developed for the coupling of the dynamics of DNA 
denaturation zones and selectively single-stranded DNA binding proteins |2Z1 12H] • This 
general scheme allows us, as limiting cases, to consider the different regimes of fast and 
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slow chaperone (un)binding mentioned above, including the explicit limit case of pure 
ratcheting. In that course, we concentrate on the effect of the chaperones, and we view 
the biopolymer as a stiff chain. Additional entropic effects due to the accessible degrees 
of freedom of the chain can be straightforwardly incorporated, as is briefly discussed. 
We also do not consider anomalous diffusion dynamics of very long translocating chains 
here [211 EDI IHH ■ We note that the master equation approach is a natural powerful 
approach to translocation dynamics (and similar random processes in systems with a 
discrete coordinate) , as was demonstrated in references jHHl • 

We finally point out three differences between our approach compared to the 
original Brownian ratchet ideas put forward in reference ^1^]. (1) In ^2 Ej the 
polymer motion through the pore was assumed to proceed through continuous diffusive 
motion until a chaperone binds; an implicit assumption in that approach is that the 
characteristic diffusive step length is much smaller than the size of a chaperone. In 
contrast, we assume that the motion proceeds through random jumps of length a with 
a rate k, where a is taken to be of the order the length of a base, base-pair or aminoacid. 
The two approaches are identical when all other length scales of the problem are much 
larger than a; however in our problem the size of the chaperones can be of the order 
of a, which in the general case (and in particular for the case of strong binding) alters 
the dynamics. The justification for our approach lies in that experimentally it is found 
that the mean translocation time of single-stranded DNA through an a-hemolysin pore 
is roughly 100 fold smaller than expected including only hydrodynamic friction in the 
pore |22j. Also, for protein translocation the translocation times are many orders of 
magnitude larger than expected from simple hydrodynamic considerations jTHj. These 
results suggest that there are strong interactions between the pore and the translocating 
polymer; adopting a similar explanatory description as in [23^ in which the translocation 
process proceeds through jumps in a saw-tooth-like pore potential, our rate constant 
k is simply the Kramer's rate associated with a jump in this potential and a is the 
distances between minima. (2) In reference IT3] the chaperones bind independently 
along the polymer. Here, we include the fact that typically chaperones are larger than 
the size of a monomer, giving rise to a "car parking effect" - a chaperone can only 
bind if the space between already bound chaperones equals or is larger than the size of 
the chaperones; we thus extend the model in ^1 1^ in order to include the fact that 
in general chaperones do not bind independently. (3) We here consider arbitrary ratios 
of the relevant (un)binding rate and the diffusive rate k - this contrasts the results in 
references [TH IT^ where only the cases of Brownian ratchet motion and fast (un)binding 
dynamics were considered. We will point out similarities and differences between our 
approach and that in references ^JE] as they arise throughout the text. 

2. Master equation for chaperone-driven translocation 

Let us consider the general translocation dynamics for a chain of length Ma, where a is 
the typical linear monomer (base, base-pair, aminoacid) size. The coordinates we choose 
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for the description of the chap erone- driven biopolymer passage are the number m of 
monomers of the chain on the trans side B (the number of monomers on the cis side A is 
then M—m) and the number n of bound chaperones, see figure^ The size of the binding 
interface of a chaperone bound to the biopolymer is Xa. No chaperones are assumed 
to be present on side A. Due to inherent thermal fluctuations the variables m and n 
are stochastic variables and the aim is to understand the temporal behaviour of these 
variables; below we use a master equation formulation for the associated probability 
distribution. This description rests on the assumption that all other variables in the 
system are fast in comparison to the rate of change of the coordinates m and n. From 
the solution of this master equation any experimental observable can be calculated. 
Here, the main quantity of interest is the mean first passage time t(^) a function 
of the chain length M and chaperone size A, as well as of two dimensionless parameters 
introduced in section 01 the binding strength k and the ratio 7 between the rate k 
associated with making a diffusive step and the rate q for chaperone unbinding. In our 
description, the size of the pore and the number of monomers present in the pore at a 
given time enters only through the effective rate k (see references 123 IHHl IHH IHE] for 
more detailed investigations of the effects of finite-sized pores). 

Denote by P{m,n,t\m' ,n') the conditional probability density of finding m 
monomers of the chain on side B (the already translocated distance in units of 
monomers) with n bound chaperones at time t, provided that at the initial time the 
system was in a state with m = m' and n = n' . With the short-hand notation 
P(m, n, t) = P{m, n, t\m', n'), the time evolution of P(m, n, t) is governed by the (2+1)- 
variable master equation P7j 

oPiin Ti t) 

^ ' ' ^ = t+(m - 1, n)P{m - 1, n, t) + t"(m + 1, n, t)P{m + 1, n, t) 

— [t^(m, n) + t~(m, n)] P(m, n, t) 

+ r^(m, n — l)P(m, n — 1) -\- r^(m, n + l)P(m, n + l,t) 

— [r^(m, n) + r~(m, n)] P(m, n, t). (1) 

The transfer coefficients represent the rates for a change in the translocation 
coordinate m and the coefficients specify the (un)binding of chaperones, i.e., a change 
in n; both transfer coefficients are explicitly defined in section |3] :j: Equation (P) is 
subject to the general initial condition P(m, n, t = 0\m',n') = Sm,m'Sn,n', where ^21,22 
denotes the Kronecker delta. Here, we assume that m' = and n' = 0: 

P{m,n,0) = 6m,oSn,o, (2) 

I For equation |^ to be completely specified we need to describe the transfer coefficients just outside 
the configuration lattice (illustrated in figure [SJ as well: 

t-{m = M + = t+(m = jX-l,n^ n'"^''(m + 1)) = 

for integer j, and 

r+(m,7i = -1) = r^{m,n = n'°''''(m) + 1) = 0. 
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i.e., that initially, the chain is fully on side A, with its head just in the pore (and 
therefore has no bound chaperones), corresponding to the typical situation in vitro and 
in vivo. 

We now consider the boundary conditions that must be imposed on the transfer 
coefficients: Once the chain has arrived completely in side B, i.e., m = M, we assume 
that it cannot find its way back into the pore. We impose, that is, the absorbing 
boundary condition § 

t-(M+l,n) = 0. (3) 

For the absorbing condition (jH} imposed here, the Mth transfer coefficient t"^(M, n) in 
general is finite, giving rise to jumps out from the configuration lattice (see Figure |21), 
such that in this case J2mn^i^^^^^) decreases with time, and therefore P{m,n,t) is 
an improper probability distribution. At m = 0, we impose the reflecting boundary 
condition 

t-(0,0) = 0. (4) 

This condition guarantees that the chain does not retract from the pore to side A. 
Additionally, we require that the chain cannot move towards side A, when it is fully 
occupied with chaperones: 

t" (m = n°'^^(m)A, n°'^^(m)) = 0, (5) 

where 

n°^"^(m) = [m/A] (6) 

is the maximum number of bound chaperones for a given m (and [■] is the Landau 
bracket returning the integer value of the argument). The boundary condition (jSJ is 
crucial for the emergence of ratchet motion as discussion in subsection 15.21 For the 
(un)binding rates, we have the natural condition 

r-(m,0) = 0, (7) 

stating that no further unbinding can occur once the chain is empty of chaperones. 
Similarly, chaperones cannot bind to a fully occupied chain: 

r+ (m,n'"^^(m)) = 0. (8) 

The configuration lattice on which the jump dynamics, as defined by equation ((T)) and 
the transfer coefficients and r^, is schematically shown in figure El together with the 
boundary conditions above. 

For internal points on the lattice in figure Elthe transfer coefficients are constrained 
(but not fully determined) by the detailed balance condition ^7\. For our specific 
system, the detailed balance condition takes on the structure (compare to references 

123 EH!) 

f^lm — 1, n)S'{m — 1, ra) = t~(m, n)^{m, n) (9) 

§ We note that a reflecting boundary condition at m = M, as occurring for instance in the description 
of DNA- breathing dynamics [SZIEHIj would have the structure t"'"(M, n) = 0. 
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Figure 2. Configuration lattice spanned by translocation coordinate m and number 
n of bound chapcroncs, with allowed configuration states of the system (□). Each 
time m reaches a multiple of A, an additional chaperone can bind. The scheme also 
indicates forbidden jumps. 



and 

r~^{m,n — l)iF(m, n — 1) = r~{m,n)^{m,n). (10) 

where 3f{m, n) is the partition coefficient for a given m and n. Whereas equation (jlUj) 
holds for all n, equation Q only holds for < m < M, but not for m = M + 1. This is 
due to the fact that t"'"(M, n) ^ corresponding to our absorbing boundary condition; 
in the case of a reflecting boundary condition, t~^{M,n) = 0, detailed balance would be 
fulffiled for all m, and for long times P{m, n, t) would reach the stationary density: 

P'\m, n) = ^{m, n) /3f (11) 

where ^ = Xlm n ■^("^' Here, in contrast, our absorbing condition provides that 
P{m, n,t oo) = 0, and P{m, n, t) is in fact an improper probability distribution. An 
explicit expression for iF(m, n) and for physically realistic transfer coefficients, satisfying 
the detailed balance conditions above, are given in section HI 

The general solution to the master equation (P) can be decomposed into the 
eigenmodes according to 

P(m, n, t) = P{m, n, t\m', n') = ^ Cp(m', n')Qp{m, n)e"''''*, (12) 

p 

where the expansion coefficients Cp{m',n') are determined by the initial condition. 
Inserting above expansion into equation ((H) produces the eigenvalue equation 

t^{m — 1, n)Qp{m — l,n) + t~(m + 1, n)Qp{m + 1, n) 

— [t^(m,n) + t~{m,n)~\ Qp{m,n) 

+ r+(m, n — l)Qp{m, n — 1) + r~(m, n + l)Qp{m, n + 1) 

- [r+(m, n) + r~(m, n)] Qp{m,n) = -ripQp{m,n) (13) 

for the pth eigenmode, with eigenvalues rjp and eigenf unctions Qp{m,n). We label the 
eigenvalues such that < 770 < ''7i < ''72 • • • < Vm- We note that all rjp are real and 



positive (see Appendix A and [SHI EH]) which guarantees that P{m,n,t — * 00) = 0, 
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as it should for an ergodic system with a sink (see equation Q). The eigenvectors 
Qp{m,n) satisfy the orthonormahty relation given in equation ()A.4|) . We prefer using 
an eigenvalue approach (spectral representation) to the present problem rather than 
solving the master equation in real time, since the eigenvalue approach avoids time 
discretization problems. 



3. Mean first passage time of translocation 

We here derive an expression for the mean first passage time of translocation, i.e., the 
mean time it takes the biopolymer, while being driven by chaperone (un)binding, to 
fully cross the membrane pore. We will determine the mean first passage time as a 
function of the eigenvalues rjp and the eigenfunctions Qp{m, n) defined through equation 

m- 

As before, denote by P(m, n, t\m', n') the conditional probability to find the system 
in state {m,n) at time t, given the initial condition {m',n') at time t = 0. Then the 
(survival) probability that the absorbing boundary at m = M + 1 has not yet been 
reached up to time t is 

o5^(m', n', t) = P{m,n,t\m' ,n'). (14) 

m=0 n=0 

The probability that the absorbing boundary is reached within the time interval [t, t + dt] 
is 

f{m', n', t)dt = ,y{m', n', t) - ,y{m', n', t + dt) 
d 



^^y{m',n',t)jdt. (15) 

This expression is positive, as y is decreasing with time. From the first passage time 
density f{m',n',t), the mean first passage time thus follows by integration: 

/»00 /"OO 

j{m',n')= / tf{m',n',t)dt= / y{m',n',t)dt. (16) 



To express the mean first passage time j in terms of the eigenvalues and eigenfunctions, 
we introduce the eigenmode expansion (fT^ into equation (fTH|). yielding 



J (m', n 



/oo 
Cp(m', n')Qp{m, n)e~'^^^dt 
,,,,,,, - (1-1 



m,n " p 



where P^^{m^ n) is given in equation (lllj). We have above made use of the orthonormahty 
relation ()A.4|) in order to express Cp{m',n') in terms of the initial probability density 
P{m,n,0\m',n'), and used the fact that this general initial condition takes the explicit 
form P{m,n,0\m',n') = 6m,m'Sn,n'- Equation (fTTj) is the discrete counterpart of the 
continuous result derived in [38j, and expresses the mean first passage time (for any given 
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initial condition, specified by m' and n') in terms of tlie eigenvalues and eigenvectors of 
equation (fT^ . 

Using the general result ()17j) together with the special initial condition 
considered here, we see that the mean translocation time for a polymer starting with 
its head in the pore becomes: 

j = j{o,o) = J2jp%\ (18) 

p 

with 

T. = $|^E«p(".«)- (19) 



m.n 



The mean translocation time is thus obtained by summing up the relaxation times rj^ ^ 
with weight-factors jp, in turn being determined by the eigenfunctions Qp{m,n). 



4. Partition function and transfer rates 



The partition coefficient 3f{m, n) appearing in the detailed balance conditions (jHI) and 
()10|) constrain, but do not fully specify, the transfer coefficients t^(m, n) and r^(m, n). 
In this section we give explicit expressions for both the partition coefficient as well for 
the transfer coefficients. 

For a given value of translocation coordinate m and a given number of bound 
chaperones n the partition coefficient is: 

J^(m,n) = J^P°^y(m)^'=^^P(m,n), (20) 

corresponding to the product of the partition coefficients S'^°^^{m) of the undressed 
biopolymer and S'^^'^^{m^n) of the chaperone (un)binding, which in turn depends on 
the translocation coordinate m defining the number of accessible binding sites, see figure 
E The polymeric factor ^P°'^(m) measures the accessible degrees of freedom (or better, 
their confinement) of the biopolymer chain threading through the pore. For a flexible 
chain, this would invoke the critical exponent for a polymer grafted to a surface |24j . 
Single-stranded DNA or RNA would quahfy for such a description. Double-stranded 
DNA and proteins, are more of a semiflexible nature, and corrections to the scaling 
behaviour would ensue. Moreover, the interactions between chain and pore are not 
fully known, compare the discussions in references fOl 01] • For simplicity, and to solely 
concentrate on the chaperone effects, we assume our chain to be completely stiff, so that 

J^P°iy(m) = 1. (21) 

Note, however, that any expression for ^P°'y(m) can easily be implemented through 
the numerical formalism discussed below, provided that the chaperone binding does not 
depend on the curvature of the polymer. 

The second factor in the partition function ^(m, n), the contribution S'^^^'^{m,n) 
from chaperone (un)binding, 

J^"'^^P(m, n) = /€"fi=^'^P(m, n) (22) 
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combines the n-fold gain of the binding energy, entering through the factor k", with the 
number of possible ways of arranging n proteins along m biopolymer monomers, 

n^'^nm^n) = - - 1)-^ = , (23) 

\ n J nl[m — An)'. 

compare references |32l I2H1 1^ • Here, we also introduced the binding strength 

= CoK'^ (24) 

that, in turn, depends on the chaperone concentration cq in the surrounding solution and 
the equilibrium binding constant K'^'^ = fgexp (|_Ebind|/^B^) with the typical volume Vq 
of the chaperones; -Ebind is the chaperone binding energy, /3 = l/ksT with ks being the 
Boltzmann constant and T the temperature. 

The transfer rates t^(m, n) and r^(m, n) governing the transitions m — > m ± 1 and 
n ^ n ± 1 in the master equation are subject to the detailed balance conditions © 
and ([rUI). Here, we seek concrete expressions for the rates. For the forward translocation 
rate t^, we solely observe the influence of the polymeric degrees of freedom. We assume 
that the rate limiting step for the polymer motion originates from interactions in the 
pore |22j and choose the form 

+ . X ,J^P°^y(m + l) , , 

t+ m, n) = k—-\——^ = = t+ m . 25 

It simplifies to the constant rate k due to the assumption of a stiff chain, j!°P°'y(m) = 1, 
where k is the rate for making a diffusive jump of length cr, see introduction and figure 
d Conversely, for the backward rate, we specify 

, ,J^(m-l,n) , ^^^^P(m- l,n) , fi^^^P(m - 1, n) , , 
t" m,n = k—— = fc— — ^ = fc— — ^ V^. 26 

In the case of the t^, we thus choose that while the forward rate is independent of 
the number of bound chaperones and their configurations, the sliding of the chain 
toward the cis side A of the membrane will, in general, be opposed by bound proteins, 
t~{m,n) < k. We have chosen t~ such that it is proportional to the probability 
^chap(^ — 1, n)/Q^^^'^{m, n) that the binding site closest to the pore on side B is vacant. 
In the limit n = 0, both forward and backward rates are identical: t"'"(m, n = 0) = 
t~(m, n = 0) = fc, and we have pure, unbiased diffusion with a bare diffusion constant 
Do = (T^/c, II as it should be. We point out that whenever the polymer is not fully 
occupied, the binding proteins may, in general, slide along the DNA thereby 
possibly exploring the different bound chaperone configurations on a faster timescale. 
Such sliding motion was present in the Brownian simulations in reference jTTj. 

For the transfer coefficients associated with chaperone (un)binding we specify: 

, , J^=i>^P(m,n + l) 
r (m, n) = (n + l)q — — -r — -, : — 

II The bare diffusion constant Dq sfiould not be confused witfi tlie effective diffusion constant relevant 
for general driven motion over distances much larger than a, see reference |23| . 
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(n + l)fi=^^P(m,n + l) 



(27) 



r (m, n) = nq = k'yn. 



(28) 



To make comparison with the chaperone unbinding rate q easier, we introduced here 
the dimensionless ratio 



between the the unbinding rate g of a single chaperone (for n = 1 we have r~ = q) 
and the diffusive rate k. We have chosen r+ to be proportional to the number of ways 
^ = {n + l)Q'^^^^{m,n + 1)/Q^^^'^{m,n) of adding one additional chaperone if there 
are already n bound. In the general case ^ includes parking effect" |ini2T|. 

i.e., the fact that a chaperone can only bind in between two chaperones which are 
separated by a distance equal to or larger than the chaperone size. In the limit A = 1 
we have independent binding and the resulting expression becomes ^ = m — n as it 
should, and r"*" is thus proportional to the number of free binding sites, compare the 
discussion in |2Hl IHZ] • We note that r"*" depends hnearly on the volume concentration cq 
of chaperones (compare reference ^1]) through the parameter k = cqK'^'^, where K^^ is 
the binding constant as before. Conversely, the unbinding rate is simply proportional to 
the number of chaperones bound to the biopolymer (as it should jHZI) and independent 
of the concentration of chaperones. A word on the choice of this asymmetric form is in 
order. To put the full effect of chaperone concentration and binding energy contained 
in ^f^^^p into the binding rate r"*", and therefore leave the unbinding rate r~ unaffected 
corresponds to the physical picture that the unbinding process should solely depend 
on the dissociation of local bonds. It should also be noted that cooperative binding of 
chaperones can in principle be included in our formalism [21^. 

The explicit expressions (|^ . pHjl . ^I7\ and (j2Hl) for the transfer coefficients, 
together with the boundary conditions in equations ©, ©, © and (jHJ completely 
determine the translocation/chaperone (un)binding dynamics through the master 
equation ((H) or the associated eigenvalue equation (fT^ . In the next section we solve 
these equations in order to obtain the mean translocation time, using the expression 
fll8|) derived in the previous section. 

5. General and limiting results 

In this section we obtain results for the mean translocation time for different physical 
parameters. In particular three cases can be distinguished: (i) slow chaperone binding 
leading to unbiased pure diffusion for a finite translocating chain; (ii) fast chaperone 
binding but slow unbinding leading to a pure chaperone-mediated ratchet motion; and 
(iii) fast chaperone (un)binding dynamics tractable via adiabatic elimination, leading 
to a one variable effective force, compare figure El This analytical treatment is 
complemented with (iv) numerical analysis. We stop to note that the emergence of 
the three limiting cases from a general scheme is the major advantage of the approach 




(29) 
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pursued herein in comparison to the one-dimensional treatment in reference which 
exphcitly accounts for case (iii) above only. 

5.1. Slow binding dynamics, 7/? — > 

Slow chaperone binding is prevalent when the dimensionless ratio 7K — 0, as in this 
limit we have that r~^{m,n) (see equation ()27|)). Experimentally one can reach 
small values of k, and hence the slow binding regime, by sufficiently low chaperone 
concentration cq or for small binding energy -Ebind? see equation ()24|) . As initially n = 
(see equation (j21)), the dynamics is fully described by the transfer coefficients 

t+(m,0) = t~(m,0) = A;, (30) 
together with the reflecting boundary (compare equation (jH)) 

t-(0,0) = (31) 
and the absorbing boundary condition (see equation (jH))) 

t-(M + l,0)=0. (32) 

We see from equation ()3()j) that in the slow binding limit the rate for a forward step 
equals the backward rate and therefore the motion is purely diffusive for a finite chain. 
Making use of the general expression for the mean first passage time for one- dimensional 
motion contained in equation (jB.lll we arrive at 

M m 

T = Tdiflf = XI XI = 

m=0 m'=0 

For large M, the mean first passage time thus scales like j ~ M^/2/c, as it should for a 
purely diffusive first passage process. In terms of the bare diffusion constant Dq = a^k 
and the chain length L = Ma we have j ~ L?'/2Dq in agreement with The 
diffusive motion on the general configuration lattice is schematically illustrated by the 
blue dashed arrows in figure El (top). Notice that due to the slow binding only a small 
part of the phase space is explored before the polymer is absorbed. It is sometimes useful 
to picture the dynamics as jump motion in a force field; for purely diffusive motion we 
have a force F{m) = 0, i.e. the force is zero, as schematically illustrated in figure El 
(bottom). In figure El the mean translocation time (j33|l is shown and compared to the 
different cases discussed in the remaining part of this section. 

5.2. Fast binding but slow unbinding, 7 — >■ 0, 7K — > 00 and A = I.' ratchet motion 

As the effective rate 7 and the binding strength can be chosen independently, we can 
specify them such that we encounter slow unbinding designated by 7 ^ but fast 
binding, 7/? — > 00. For small chaperones, A = 1 (univalent binding), this corresponds to 
the stepwise transitions marked by the red arrows in figure El Each time a monomer of 
the biopolymer exits from the membrane pore towards the trans side, it is immediately 
occupied by a chaperone. The irreversibly bound chaperones do not unbind and prevent 



(M + 1) 



(33) 
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Figure 3. Top: Schematic of the analytically tractable limiting cases of the coupled 
translocation-chaperone (un)binding dynamics, (i) The dashed blue arrows indicate 
the purely, unbiased diffusive motion for slow binding, (ii) The solid red arrows indicate 
the translocation dynamics for ratchet motion in the case of fast chaperone binding. 
Notice that for these two cases only a small part of the phase space is explored before 
the polymer is absorbed at m = M + 1. In contrast, for (iii) the fast (un)binding case, 
the entire phase space is explored before absorption. Bottom: Effective force F{m) 
for the three limiting regimes above. For (i) diffusive motion we have F{m) = (blue 
dashed line); for (ii) ratchet motion we effectively have F{m) ~ 2Fo (red solid line 
referring to the right ordinate), with a characteristic force Fq = ksT/a; for (iii) fast 
(un)binding we have F{m) — 2Fo{t'^{m) — t^(m))/(t"'"(m) +t~(m)) where t^(m) are 
defined in subsection 15.31 We here use a binding strength k — 0.8, and two different 
chaperone sizes A = 1 (the upper dash-dotted line) and A 5 (the lower dashed curve) . 
Note that the onset of a non-zero force is at m = A, i.e. the force is zero unless there 
are sufficiently many monomers to accommodate at least one chaperone. Also notice 
the oscillations in the force for the case A = 5. We point out that the forces above 
are not explicitly used in the dynamical scheme, but are convenient for illustrative 
purposes. 
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backwards motion (see equation so that the motion becomes unidirectional with 
effective rates V~{m) = k and t~(m) = 0. The general master equation (Q) then reduces 
to the effective (l+l)-variable equation 

= k[P{m- 1, t) - P{m, t)] , (34) 

where we introduced the short-hand notation P{m,t) = P{m,n'^^^{m),t). The process 
is that of taking M + 1 forward steps with a rate per step k, and the mean first passage 
time therefore becomes 

T = Tratchct = A;~'(M + 1). (35) 

Thus, J scales linearly with M for large M in this ratchet limit. The mean translocation 
time (jH^j) is shown in figure Equations of the same form as (jMj) also occur in 
the theoretical description of shot noise, and correspond to the forward mode of the 
wave equation. For ratchet motion we define an effective force F{m) = 2Fo(t^(m) — 
t~(m))/(t'''(m) +t~(m)) = 2Fo, where Fq = ksT/a; the effective force associated with 
ratchet motion is thus constant, as schematically shown in figure 01 (bottom). 

Comparing equations and we see that jratchet /idifr = 2/M for large M. 
This results contrasts with the results of a Brownian ratchet ^| |^ for which the same 
quantity equals 1/M. This difference originates from the fact that here the size of a 
chaperone equals the diffusive step length a, whereas in the Brownian ratchet a —>■ 
(continuous diffusion versus our stepwise diffusion), but with a finite chaperone size. ^ 

5.3. Fast (un)binding dynamics, 7 ^ oo and 'yn 00 

Finally, we address the case of fast binding, 7/? 00, and fast unbinding, 7 ^ 00. Under 
these conditions, it is possible to adiabatically eliminate the fast variable n [23 EHl, 
resulting in the effective forward and backward rates 

t-(m)= ^Hm,n)^^. (36) 

For a stiff polymer, implying t"'"(m, n) = k (see equation ()25|l ). we obtain the following 
simple expression for the forward rate: t^(m) = k. When calculating the backward rate 
t (m) one should remember to incorporate the boundary rates as given in equation 0. 
Here, in general, ^(m) = ^^"^^ (m) ^'^^'^^ (m) , and 

3f'^-P{m)= ^ J^=^^"P(m,n). (37) 

n=0 

% The results from jl4[ [T^ can be recovered if we imagine the chaperones in the ratchet case to be fast 
binding but immobile after binding. We would then expect the mean first passage time for Af ^ A ^ 1 
to be the sum of M/X diffusive mean first passage times fc~^A^/2. Thus the ratio with the purely 
diffusive mean first passage time in equation (|33|) would become Tratchct /jdiff = A/M in agreement 
with [TllTHl . 
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Making use of the fact that the transfer coefficients satisfy the detailed balance condition: 
t'''(m — l)^(m — 1) = t~(m)iF(m) (see reference j2H]) so that we can use equation ()B.3|) . 
we obtain the mean first passage time 

' ^,t+(m)^(m) ^n^M4l, 

m=0 ^ ' ^ ' m'=0 m=0 m'=0 

where we have used that here t^(m) = k. The sum in (jH^ can straightforwardly be 
numerically evaluated and the resulting mean translocation time is shown in figure 

In order to simplify the mean translocation time ()38j) further for large chain lengths 
M, a knowledge of the partition function ^{m) for large m suffices. Instead of the 
combinatorial approach given in the previous subsections, the partition function for 
large m is obtained using the approach in reference 01j (see also jlSj), giving 

^ CA- (39) 

where C is a constant and (assuming no cooperativity effects IIUEII) Amax is the largest 
root A to the algebraic equation 

A^ - A^-^ - K = 0. (40) 

The equation above is of order A, i.e., the order is determined by the size of the 
chaperones. Introducing equation (jH^ in equation (PH| we find that the mean first 
passage time for very large M approaches 

T |a/»i - Tad = k-' (M + 1), (41) 

-'^max ^ 

which proves that for large chain lengths the adiabatic result for j scales as M , just as 
for ratchet dynamics, but in general with a smaller prefactor than for ratchet motion. 
For the case of weak binding, n small, so that InA^ax ^ 1 we can write equation ()4ip 
according to j ~ k~^ In A~^^(M + 1), so that j oc In A~?^^ in agreement with the results 
in reference pTj . 

For univalent binding (A = 1) we solve equation (jlUj) and find the explicit result 
Amax = (1 + k). Interestingly, the mean translocation time is then inversely proportional, 
J oc 1//, to the filling fraction, f = k/{1 + n), of the chain, see appendix C in reference 
PT] : since < / < 1 the translocation process is slower than for ratchet motion as 
it should be, compare equation (jH^jl . We note that Tad/ldifr = 2(1 + k~^)/M, see 
equations (jH^ and (jH}. + For divalent binding (A = 2) equation (pOj) becomes 
a second order algebraic equation, which can straightforwardly be solved, yielding 
Amax = [1 + (1 + 4k)^/^]/2. In figure El the result as contained in is shown for 
the cases A = 1 and A = 5; for the case A = 5 the algebraic equation is solved 
numerically. 

We point out that, in contrast to the cases of slow binding and ratchet motion 
discussed in the previous subsections, for the fast binding and unbinding case (reversible 

+ The corresponding result for a Brownian ratchet is Tad/ldiff = (1 + 2k^^)/M . The difference has the 
same origin as discussed in subsection 15.21 
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Figure 4. Enumeration scheme for the numerical analysis: The two-dimensional grid 
points (m, n) are replaced by a one-dimensional running variable s. See text for details. 



binding [21]) considered here typically the larger part of the phase space is explored, 
see figure El (top), before absorption. For fast (un)binding we define the effective force 
F{m) = 2Fo{t~^{m) — t~(m))/ (t+(m) + t^{m)), with Fq = ksT /a as before. The force is 
illustrated in figure El (bottom), showing a constant value for large m; this asymptotic 
value is smaller than that for ratchet motion. Also, notice that for A = 5 the force 
oscillates with a period ~ A for smaller m in agreement with reference [21]. 

5.4- General case, numerical analysis 

To solve the eigenvalue equation (fT!?|) by a numerical scheme, it is convenient to replace 
the two-dimensional grid points (m, n) by a one-dimensional coordinate s counting all 
lattice points, compare and jlOI. We choose the enumeration illustrated in figure IH 
From this figure we notice that n e [0, A^] (where = [M/A]) and that m G [nA,M]. 
An arbitrary s-point can be obtained from a specific (m, n) according to: 

1) 



s = n{M + 1) -^n'X + m = n{M + 1) 



n(n 



-A -|- m. 



(42) 



n'=l 



From this relation we notice that the maximum s value is 
S = maxjsl = (AT + 1) ( M + 1 - ^ 1 - 1. 



max{s} = {N + 1) (^M +1 - -1. (43) 

Expression ()42|1 allows us to change the transfer coefficients to the s- variable, t^(m, n) —> 



t^(s) and r^fm, n) 



[s), using the explicit expressions (j^ . (j^ . (j2I|) (^Hj) for the 



transfer coefficients, together with the boundary conditions in equations (jlj), ©, 
and (|H1). From equation and figure 01 we notice that a local jump in the m-direction 
corresponds to a local jump also in s-space, i.e., that 



m ^ m-1 s ^ s - l|i<m<M , 
m ^ m + 1 <^=^ s ^ s + l|o<m<M-i • 
However a jump in the ra-direction is equal to a non-local jump in s-space: 

n — ^ Tl — 1 \ ] ' S — ^ S — Z\s I l<n<n"^^^ (m) 5 
n ^ n + 1 <^==^ S ^ S + As^\o<n<n^'^^{m)-l , 



(44) 



(45) 
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with As = M + 1 — nX and As"*" = M + 1 — (n + 1)A. Using the results above, we find 
that the eigenvalue problem (fT^ can be written in matrix form as 

J2w{s,s')Qp{s') = -rjMs) (46) 

s' 

where explicitly the matrix-elements are 

W{s,s-1) = t+(s - 1), for s dl 1 < m < M, 

= 0, for s iti m = riA A n = n^'^'^{m) 
W{s,s + l) =t-(s + l), forsrhO<m<M-l 
W{s, s - As") = r+(s - As"), for s ^ 1 < n < n^^^(m), 
W{s, s + As+) = r+(s + As+), for s rh < n < n"^^"(m) - 1, 
W{s,s) = -(t+(s)+t-(s) + r+(s) + r(s)), (47) 

and the remaining matrix elements are equal to zero. We have introduced the notation 
s rh with the meaning "s is to be taken for" . Notice that (with the present enumeration 
scheme) the tridiagonal part of the l^-matrix is determined by jumps in the m-direction, 
i.e., by t+ and t~, whereas elements outside of the tridiagonal part are determined by 
r+ and . The problem at hand is that of determining the eigenvalues and eigenvectors 
of the (5 + 1) X (5* + l)-matrix W above. Convenient checks of the numerical results 
include: (i) The eigenvalues should be real and negative (so that rjp > 0); (ii) The 
eigenvectors should satisfy the orthonormality relation, equation ()A.3|) . 

In figure El we show results for the mean translocation time. From the results of 
the previous section we see that for ratchet motion (subsection 15. 2|) as well as for fast 
binding and unbinding dynamics (subsection I5.3|l the mean translocation time scales as 
~ M. For slow binding as discussed in subsection 15.11 we have that j scales as ~ 
for large M. The general results shown in figure El are obtained by numerically solving 
the eigenvalue equation ()4fjj] . and using the mean first passage time expression (fTTH) . All 
results are shown for for the case of small chaperones A = 1, and for larger chaperones, 
A = 5. We notice the slow binding result (jHSl) and the ratchet results constitute 
upper and lower limits to the general mean translocation time j. Also, the approach to 
the very large M result PT|) is rather slow, due to finite size effects. 

We point out that even though our main focus here has been to calculate the 
mean first passage time j, any quantity of interest (such as, for instance, the mean 
first passage time density fit), see section 0} can be obtained from the solution of the 
eigenvalue problem presented in this section. 



6. Conclusions and discussion 



In this study the translocation of a stiff polymer consisting of M monomers through 
a nanopore in a membrane was investigated in the presence of binding particles 
(chaperones) that bind onto the polymer. Our (2+l)-variable master equation governing 
the coupling of the random biopolymer passage and the random chaperone (un)binding is 
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Figure 5. Mean translocation time j as a function of AI, using (top) linear axes; 
(bottom) logarithmic axes. Results are shown for: (i) slow binding, equation (ii) 
the ratchet result, equation (|35|l : (iii) fast (un)binding, equation H38|) for different sizes 
of the chaperones (here A = 1 and A = 5 respectively) with binding strength k = 0.8. 
Also, results are shown for (iv) the general dynamics, as obtained by numerically 
inverting the transfer matrix (|47|l and using the mean first passage time expression 
(|19|l . In addition, in the bottom graph the large M result H41|) is shown. 

an extension of a previous continuum Fokker-Planck approach JT] . From our approach 
three hmiting regimes emerge: (i) For the case of slow binding the motion is purely 
diffusive, and the mean translocation time j scales as j ~ for large M; (ii) for 
fast binding but slow unbinding, the motion is ratchet-like, for small chaperones A = 1, 
and J ~ M; (iii) for the case of fast binding and unbinding dynamics, we performed 
the adiabatic elimination of the fast variable n, and find that for a very long polymer 
J ^ M, but with a smaller prefactor than for ratchet-like dynamics. We point out that 
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our approach is the first, to our knowledge, that is able to account for (the approach 
towards) these three limiting regimes within a single framework. For the intermediate 
cases our scheme is solved numerically, which requires numerically inverting the transfer- 
matrix W (of approximate size x M^), for which an explicit expression is given. 

Let us summarize the essential differences between our approach compared to the 
original Brownian ratchet ideas in reference ^1^]. (1) In this reference the polymer 
motion through the pore was characterized by continuous diffusion. In contrast, we 
assume that the motion proceeds through random jumps of finite step length, as 
appropriate, for instance, for a polymer moving through a saw-tooth like pore potential, 
as suggested in |^ for explaining the surprisingly large translocation times for DNA 
translocation. (2) In reference the chaperones bind independently along the 

polymer. Here, we include the fact typically chaperones are larger than the size of a 
monomer, giving rise to a "car parking effect". (3) We here consider arbitrary ratios of 
the relevant (un)binding rate and the diffusive rate - this is in contrast to the results in 
references |14[[TK] where only the cases of Brownian ratchet motion and fast (un)binding 
dynamics were considered. 

An implicit assumption in our study is that the number of binding proteins in 
the bath surrounding the translocating chain is large (but of finite concentration), 
see reference [25, appendix B. In other words, we assume that the diffusion of 
the chaperones in volume is fast in comparison to the time scales of translocation 
and chaperone (un)binding dynamics (well-mixing condition). Thus, chaperone bath 
depletion effects, that may occur for a finite number of chaperones ^I], are not present 
in this study. We note that we have not explicitly included chain flexibility into the 
present scheme. We have also assumed that the binding energy is the same along the 
polymer, and that the major friction for the polymer originates from interactions in the 
pore (we hence neglected hydrodynamic friction from the polymer part sticking out of 
the pore). Although each of these effects may be relevant to experimental situations 
of chaperone-driven translocation, here we focused on the pure chaperone effect. Note 
that one efficient way to include additional details in the translocation scenario is the 
stochastic simulation technique based on the Gillespie algorithm [1^1 HZ]- It would 
also be interesting to compare the results from this study with Brownian dynamics 
simulations along similar lines as in reference ^T], but with the (realistically) large 
chaperones considered here. 

We finally point out that the coupled dynamics studied here is a generic example of 
how a stochastic process can rectify another one. Our (2+l)-variable master equation 
formalism provides a general framework for treating such rectification processes, that 
may also be useful in other fields. 

Appendix A. Proof of the orthogonality and the negativeness of eigenvalues 

In this appendix we prove an orthonormality relation for the eigenvectors Qp{s) and 
that the corresponding eigenvalues —rjp are negative, see equation 
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In terms of the running variable s, see equation (|42j) . and the TV- matrix defined in 
equation (jTTj) the detailed balance conditions Q and (fTUj) can be written: 

W{s, s')P'\s') = W{s\ s)P'\s), (A.l) 

for s 7^ s'; we note that the equation above holds trivially also for s = s'. Here, 
P^*(s) = ^{s)/^ with J^(s) being the partition coefficient and ^ = J2s ^{s). 

To derive the orthogonality relation from the detailed balance condition (jA.lll . 
consider the expression 

K'lp 'Ip') P^^(s) 

s ^ ' 

^ sr^ Wjs, s')Qpis')Qpis) ^ Qp{s)W{s, s')Qpis') 

s,s' s,s' 

^ Wjs, s')Q,is')Q,is) ^ Wjs', s)Q,is)QAs') 

s,s' s,s' 

where we used equation (jlHl), and the detailed balance condition ()A.1|) . Therefore, the 
orthogonality relation 

^ Qp{s)Qp'{s) 

1^ — = ^.p'' 

s 

follows. Returning to the double index {m,n), this is equivalent to the statement 

^ P^^) = ^^-^^ 

m-,n ^ ' 

for the eigenfunctions Qp{m,n). 

To prove that all the eigenvalues are negative it is convenient to separate the W- 
matrix into the corresponding refiective matrix W^{s,s') and the additional absorbing 
terms 

W{s, s') = W,{s, s') - w,{s)6s,s', (A.5) 
where the absorbing part is 

s' 

It follows from equation (jTfj) that 

w.is) > 0, (A.7) 
with at least one of the Wa{s) strictly greater than zero, and 

Y,W,{s,s') = 0, (A.8) 

s 

which is the condition that probability is conserved when there are no absorbing terms. 
Note that it does not matter for the detailed balance conditions, equation (jA.lll . that 
the diagonal of the IV-matrix is modified. Therefore detailed balance is also satisfied for 
the refiective matrix: Wr:{s, s')P^^{s') = Wt:{s', s)P^^{s). This detailed balance condition 
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together with equation ()A.8|) guarantees that P^^{s) is an eigenvector of the M/r-matrix 
with zero eigenvalue ^7\. Furthermore, since the system is ergodic, i.e. it is possible 
for the system to reach any state from any other state, any eigenvector of VFr with zero 
eigenvalue will be proportional to P^^{s) jSZ]. The negativeness of the eigenvalues —rjp 
follows then because 

^ Wis,s')Q,is')Q,is) 



Vp 



s,s' 



_ Wrjs, s')Qp{s')Qp{s) ^ ^ n (AQ\ 

The first term after the last equality sign is non-positive since the VFr-matrix is negative 
semi-definite jHZl- The only possibility for this term to be zero is that Qp{s) is 
proportional to P^^{s). In this case, however, the last term will be strictly negative 
because at least one of the wJs) is non-zero. 



Appendix B. Mean first passage time 

For a one-variable master equation with reflecting boundary condition at m = and 
absorbing boundary condition at m = M+1, and with transfer coefficient t^{m) {t~{m)) 
in the forward (backward) direction, the mean first passage time is given by the explicit 
expression (see reference jSH], chapter 7.4) 

M m 

T = V $(m) V , , , , , (B.l) 



m=0 m'=0 



where 



*('») = ft 1t7§- (B-2) 

j=l VJ J 

and we have assumed that initially m = 0. If we invoke the detailed balance condition 
t+(m — l)^(m — 1) = t^(m)^(m), for m = 1,...M, where S'{m) is the partition 
coefficient, equation ()B.1|) reduces to 

M m 
m=0 V ^ V ^ m'=0 

The expressions above are used in section El 
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